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A lattice model of 3 He - 4 He mixtures which takes into account the continuous rotational symmetry 
' 0(2) of the superfluid degrees of freedom of 4 He is studied in the molecular-field approximation 

and by Monte Carlo simulations in three dimensions. In contrast to its two-dimensional version, 
for reasonable values of the interaction parameters the resulting phase diagram resembles that 
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observed experimentally for 3 He - 4 He mixtures, for which phase separation occurs as a consequence 
of the superfluid transition. The corresponding continuum Ginzburg-Landau model with two order 
parameters describing 3 He- 4 He mixtures near tricriticality is derived from the considered lattice 
^ ■ model. All coupling constants appearing in the continuum model are explicitly expressed in terms 

of the mean concentration of 4 He, the temperature, and the microscopic interaction parameters 

< 

characterizing the lattice system. 
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Spatial confinements of systems undergoing continuous phase transitions perturb the fluctuation spectrum of the 
corresponding ordering degrees of freedom. This leads to a dependence of the free energy of the systems on the 
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' distance between the confining walls which can be expressed in terms of universal scaling function. Their gradient 
, renders the so-called critical Casimir forces which are the analogues of the well known electro-magnetic Casimir or 
£^ 1 dispersion forces. 

Recent developments in the theory of Casimir forces in critical and correlated fluids [1] have provided a strong mo- 
tivation for testing them experimentally in various systems. Capacity studies of 4 He wetting films near the superfluid 
transition T\ [2] have confirmed the existence of the critical Casimir effect and for temperatures T > T\ quantitative 
agreement with corresponding theoretical predictions has been found [3-5] . Similar effects have been observed for wet- 
ting films of binary liquid mixtures near the critical end-point of their demixing transition [6] . Additional evidence for 
the critical Casimir force and detailed data have been reported for wetting films at solid substrates of 3 He- 4 He mixtures 
near their tricritical point [7,8]. These latter experiments have also raised new interesting challenges for the theory, 
which have motivated the present work. Among them is the sign and the amplitude of the Casimir force or, more 
CD ■ generally, the form of its scaling function for different values of the concentration of 3 He atoms. Theory predicts that 
these features of the Casimir force crucially depend on the type of boundary conditions which the confining surfaces 
impose on the order parameter [1]. For example, the force should be attractive for symmetric boundary conditions 
and repulsive for nonsymmetric boundary conditions. The distinction between the surface universality classes is also 
relevant. In the case of pure 4 He films near the A-point the boundary conditions at the two confining interfaces of the 
wetting layer seem to be very well approximated by symmetric Dirichlet boundary conditions forming the so-called 
ordinary surface universality class, because the quantum-mechanical wave function that describes the superfluid state 
vanishes at both interfaces [1,2]. For films of 3 He- 4 He mixture the situation is less clear. In these systems a 4 He-rich 
layer forms near the substrate-fluid interface, which may become superfluid already above the A-line [9-11], whereas 
there is an enrichment of 3 He near the opposing fluid-vapor interface. Thus the two interfaces impose a nontrivial 
concentration profile, which in turn couples to the superfluid order parameter. The experiment of Ref. [7] reports 
a repulsive Casimir force at the tricritical point but it is not immediately obvious that the concentration profile 
induces effectively nonsymmetric boundary conditions for the superfluid order parameter, i.e., symmetry-breaking 
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boundary conditions at the substrate-fluid interface (also known as the so-called extraordinary or normal universality 
class) and Dirichlet boundary conditions at the fluid-vapor interface. The superfiuid order parameter possesses a 
continuous 0(2) symmetry so that if the layer is effectively two-dimensional it is in the Kosterlitz-Thouless phase 
with the superfiuid order-parameter equal to zero [12]. On the other hand, upon approaching the A-line from the 
high temperature side the superfiuid layer thickens due to the increase of the correlation length and a dimensional 
crossover to a three-dimensional superfiuid phase with non-zero order-parameter should take place [13]. In order to 
be able to interpret and to understand the features of the Casimir force and other surface and finite-size effects in 
3 He- 4 He mixture films near the tricritical point systematic studies of a model system are needed. Due to the universal 
character of the critical Casimir force it is sufficient to choose as a model system one which belongs to the same 
universality class as the actual physical system. The prerequisite of such future studies is a detailed analysis of the 
bulk properties of such a suitable model. This is the purpose of the present paper. The model should resemble the 
main features of the bulk phase diagram of 3 He- 4 He mixtures in three dimensions (d = 3) and take into account the 
continuous rotational 0(2) symmetry of the superfiuid degrees of freedom of 4 Hc. 

The general features of phase separation and superfluidity in three-dimensional mixtures of liquid 3 He and 4 Hc 
are well known from experiments [14]. In pure 4 He there is a transition from a normal fluid to a superfiuid phase 
characterized by a complex order parameter. If 4 He is diluted with 3 He, the superfiuid transition temperature is 
depressed. Simultaneously, the tendency toward phase separation increases and at a critical 3 He concentration the 
mixture undergoes a first-order phase transition into a 4 He-rich and a 3 He-rich phase, of which only the 4 He-rich 
phase is superfiuid. In the temperature - 3 He concentration plane (T, x) the line of second-order superfiuid transitions 
T a (x) meets the boundary of the two-phase coexistence region at the tricritical point (T t w 87mK,x t « 0.67); 
x = N 3 /(N 3 + N A ), where N i7 i = 3, 4, denotes the number of atoms of 3 He or 4 He, respectively. 

In this paper we consider a simple lattice model known in the literature as the Vectoralized Blume-Emery-Griffiths 
model (VBEG) [15,16]. It is defined in Sec. II. The bulk phase diagram of the VBEG model was investigated 
only in spatial dimensions d = 2 by means of Migdal-Kadanoff recursion relations [15,16]; no tricritical point was 
found for any value of the model parameters. Here we study the three-dimensional version of this model within the 
molecular-field approximation and by Monte Carlo simulations, and we demonstrate that for reasonable values of 
interaction parameters the resulting phase diagram resembles topologically that observed experimentally for three 
dimensional mixtures, for which the phase separation appears as a consequence of the superfiuid ordering. This is 
carried out in Sec. Ill and Sec. IV, respectively. In Sec. VI we derive a two-parameter continuous Landau-Ginzburg 
(GL) model describing bulk 3 He- 4 He mixtures near a tricritical point starting from the modified VBEG model. The 
GL approach has many advantages and it is worthwhile to have a LG model with coupling constants explicitely related 
to the measured quantities, such like like temperature, superfiuid density, concentration of 3 He atoms, or parameters 
describing interactions. We close our paper with a summary and conclusions. 



We consider a simple lattice model of 3 He- He mixtures which takes into account the continuous rotational 0(2) 
symmetry of the superfiuid degrees of freedom of 4 He. It is a vectorial generalization of the spin-1 model used by 
Blumc, Emery, and Griffiths (BEG) [17] for describing the A-line and the tricritical point in 3 He- 4 He mixtures. In 
this model each simple cubic lattice site i is associated with an occupation variable ti, taking the values or 1, and a 
phase 9i (0 < 0, < 2n) which mimics the phase of the 4 He wave function. A 3 He atom at site i corresponds to U = 
and a 4 He atom to ti = 1. Since the model in this reduced version does not allow for unoccupied sites, the model 
does not exhibit a vapor phase. (In future studies the model can be generalized to incorporate the vapor phase; here 
it is left out for reasons of simplicity.) 9i reflects the superfiuid degrees of freedom. The model Hamiltonian consists 
of a lattice gas part describing a binary mixture and a term responsible for the "superfiuid" ordering. Since only 4 He 
atoms couple to the superfiuid order parameter the Hamiltonian is taken as 



where the first two sums are over nearest-neighbor pairs < ij >, and the last sum is over all lattice sites. The lattice 
constant a is taken to be equal to 1. 

In the lattice gas model of the 3 He- 4 He binary mixture the coupling constant K and the field A are related to the 
effective "He-^He interactions -K a fj [18] 



II. THE MODEL. 
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K = K 33 + K 44 - 2K 3Al 



(2.2) 



and to the chemical potentials /i 3 and \x A of 3 He and 4 He, repectively, 
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(2.3) 



where z is the coordination number of the lattice [z = 2d, where d is the spatial dimension of the system; z = 6 in 
the present case). 

In the liquid the effective interactions K a p are different for different a and (3 due to the differences of mass and 
statistics between 3 He and 4 He atoms. The coupling constant J is related [16,21] to a bare, superfluid density po(T) 
by 

J = h 2 p {T)a {d -^/m 2 (2.4) 

where m is the mass of a 4 He atom, a is the mean interparticle spacing (the lattice parameter in the lattice model). 
The superfluid density can be measured from the velocity of third sound and from the response of a torsional oscillator 
[19]; it has units of mass per unit volume (or area in two dimensions). Here we are concerned only with the case J > 
and K > 0. 

When all occupation numbers ti are equal to I, up to constants the first term in Eq.(2.1) corresponds to the classical 
XY model (the planar rotator model) for pure 4 He. Therefore, in the limit of A — > — oo the partition function of the 
model reduces to that of XY model up to a factor e KzN where N is the number of lattice sites. 

The model as defined above is known in the literature as the Vectoralizcd Blumc-Emcry-Griffiths model (VBEG). 
It was first proposed by Berker and Nelson [16] and, independently, by Cardy and Scalpino [15] to describe thin films 
of 3 He- 4 He mixtures, for which the mechanism of the superfluid transition is different from that in d = 3; there is no 
spontaneous breaking of the continuous symmetry in d — 2 [20], i.e., the order parameter does not become nonzero 
below the transition temperature. The superfluid transition in films of 3 He- 4 Hc mixtures is of the Kosterlitz-Thouless 
type [12]. 

In d = 2 the phase diagram of the VBEG model was obtained by means of the Migdal-Kadanoff renormalization- 
group method [15,16]. Its features are qualitatively similar to those observed experimentally for the corresponding 
three-dimensional mixtures, except that there is no true tricritical point for any value of the model parameters. 
The line of the superfluid transitions (A-line) is connected to the phase-separation curve by a critical endpoint at a 
temperature distinctly lower than the phase-separation critical temperature. Thus upon lowering the temperature 
the system first phase separates into two normal fluids with different concentrations of 3 He. At lower temperatures, 
there is phase separation into a superfluid phase with a low 3 He concentration and into a normal fluid with a high 
3 He concentration. 

In this paper we are interested in the corresponding three-dimensional systems. We determine the phase diagram 
of the VBEG model within mean-field approximation and by Monte Carlo simulations. 



III. MOLECULAR FIELD APPROXIMATION 



A. Free energy 

In this section we determine the phase diagram of the VBEG model whithin the molecular-field approximation. It 
is derived from the variational method based upon approximating the total equilibrium density matrix by a product 
of local site density matrices pi [21]. 

The variation theorem for the free energy reads: 

F<F p =Tr(pH) + (l/l3)Tr(p\np) (3.1) 

where F is the exact free energy and F p is an approximate free energy associated with the density matrix p; f3 = 1 /kgT. 
The minimum of F p with respect to the variation of p subject to the constraint Trp — 1 is attained for the equilibrium 
density matrix, p — e~^ n /Tr(e~^ n ). 

Whithin mean-field theory the density matrix is approximated by 

AT 

P = Po = l[pi (3.2) 

i=l 

where in homogeneous bulk systems the local density matrix pi is independent of the site i. For the Hamiltonian 
given by Eq.(2.1) the variational mean-field free energy per site is 
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— = N = -Y {Tr{UPl)) 2 

+ ATr(Upi) + (1/ (3)Tr(pi In Pi ) (3.3) 



"" - " - -(Tr(U Pl )f- - [(Tr(Ucose iPl )f + {Tr{t l s\ne i p l )f] 



where K = zK and J = zJ. To determine variational minima to Eq.(3.3) we treat the local site density pi as a 
variational function, and the best functional form in terms of ti and 6i is obtain by minimizing F po with respect to 
Pi. In this procedure, however, the connection between F Po and the Hemholtz free energy functional of the order 
parameters is not straightforward [22]. Minimizing $/N + rjTr{pi) with respect to pi and with r\ as a Lagrange 
multiplier one obtains 

Pl = e~ 0hi /Tr{e~ 0hi ) (3.4) 
where hi is the single-site molecular Hamiltonian given by 

hi = -K(Tr(tiPi))ti - J [(Tr(U cos 6i))ti cos 8i + (Tr(U sin Oi))U sin 9 t ] + At, . (3.5) 
We define the following order parameters: 

Q = 1 - x =< ti > (3.6) 

and 

M x =< ti cos 9i > M y =< ti sin 9 t > . (3.7) 

Q corresponds to the concentration of 4 He, x to the concentration of 3 He, and M x , M y are the components of the two- 
dimensional supcrfluid order parameter M = (Mi,M 2 ) with M = M% + My . Within this approximation Q(A,T) 
and M(A,T) are given by two coupled self-consistent equations: 

= I (J3JM) 

e 0(-KQ+A) +/ ( j gjM) ^ ' ' 



and 



M= f . h ^ . (3.9) 

where Iq{z) and I\{z) arc modified Bessel functions [23]. 
The equilibrium free energy $(A, T) is given by 

$(A, T)/N = |((Q(A, T)f + J -{M{A, T)f + (1//3) ln(l - Q(A, T)). (3.10) 

Most parts of the phase diagram can only be determined by solving the equations for Q and M numerically. Some 
regions, however, can be studied analytically. 



B. A-line and tricritical point 

In order to find the line of critical points on which second-order transitions from the normal (M = 0) to the the 
superfluid (M ^ 0) state take place, one needs the thermodynamic potential in terms of the order parameter M, 

A(M, A, T) = $ - MH , (3.11) 

where H is a field conjugate to M: 

The conditions for the critical points are: 
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dH d 2 H d 3 H 



<9M dM 2 ' 5M 3 



and the tricritical point is determined by 



To find if as a function of M we use the analogue of Eq.(3.9) for H ^ 0: 

M= _hjf>JM + fiB) 

e /3(-XQ+A) + 7 (/3JM + 

Since this equation cannot be inverted explicitly, we expand it around H = keeping only terms linear in H , and find 

= Ji^JM) - MIq(PJM) + Me«-^+ A ) ^ 
MI 1 (0JM)-{l/2){I o (0JM)-I 2 (/3JM)y 

Applying the conditions formulated in Eq.(3.13)-(3.16) yields the whole line of critical points, i.e., the A-line T s (x): 

T M =l^ = f. (3 , 7) 

It follows that, as the concentration of 3 Hc increases from zero, T s decreases linearly from T„(0) = J/2. The critical 
curve A = A S (T) in the (A, T) plane can be obtained by first solving Eq.(3.8) for A (here and in the following we 
include ks into T) which gives 

A(Q,T) = Tln(l -Q) - T\nQ + KQ + TlnI (JM), (3.18) 

and then evaluate Eq.(3.18) for M = and Q = 2T/J (see Eqs.(3.6) and (3.17)). 

The line of second-order phase transitions ends at the tricritical point (T t ,x t ) where the transition changes to a 
first-order one. From Eqs. (3.14) and (3.16) one obtains for the temperature T t 

and for the concentration x t 

T t /T s (0) = 1 - x t , (3.20) 

provided d 5 H/dM 5 > holds for the chosen value of K/J. It is possible that not all critical points on this so-called 
critical line represent equilibrium phase transitions because the latter ones are preempted by first-order demixing 
transitions. Thus it can be that only a portion of this so-called critical line gives the A-line, the rest being metastable. 



C. Demixing 

For the disordered phase with M — one can easily find the first-order phase separation line from the 3 He-rich 
'normal' fluid to the 4 He-rich 'normal' fluid. The phase separation is associated with an instability loop including a 
range of Q values for which dA/dQ > and the critical point is given by dA/dQ — d 2 A/dQ 2 = 0. These last two 
relations together with Eq.(3.18), evaluated at M = 0, are satisfied if Q c = 1/2 and T c = K/A. The critical value A c 
of A is K/2. Inserting A = A c and M = into Eq.(3.18) gives 

rini— --(l/2)/?(2Q-l) = (3.21) 

For T < K/A this equation has pairs of solution (Q, 1 — Q). For M = 0, i.e., above the critical line in the (Q,T) 
plane, these solution form the coexistence curve which is symmetric about Q = 1/2 or x = 0.5. For temperatures 
lower than the intersection temperature Tj of the critical line with the curve given by Eq.(3.21) the phase rich in 4 He 
becomes superfluid and Eq.(3.21) no longer represents the coexistence curve. 
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In order to find what types of phase diagrams the present model provides we look for the phase-separation instability 
on the critical curve as determined in Subsec. Ill B and how it is located with respect to the intersection point 
Pi = (Qii Ti). Depending on the ratio K/J there are three possibilities which give three different types of the phase 
diagram: (i) the instability point P t lies below the intersection point Pi, (ii) P t lies above Pi, and (iii) the critical 
point of the transition between 3 He- and 4 He-rich 'normal' fluids falls into the instability range initiated at Pt- 



A sufficient condition for an instability loop leading to phase separation is dA/dQ > 0. Using Eq.(3.18) and the 
relation Q = MIq (f3JM) / h (j3JM) one finds 



where Q* is the critical value of Q for supcrfiuid ordering given by Eq. (3.17). Thus (<9A/<9Q)q = q„_ > 
(dA/dQ)g = Q, + and the instability will occur on the ordered side of the critical curve when (<9A/<9Q)q = q„ + = 0. 
From Eq. (3.22) the coordinates of the instability point P t are given by Eq. (3.19), i.e., they are exactly the same as 
as those of the tricritical point. We find numerically that Pi and P t = (Q t ,T t ) coincide for K/J « 2.01681. 

For K/J > 2.01681 case (i) is realized, i.e., the instability point Pt lies inside the coexistence curve between 
the 3 He-rich 'normal' fluid and the He-rich 'normal' fluid. We have obtained numerically the phase diagram for 
K/J = 2.4. It is shown in Figs. 1(a) and 1(b) in the (A, T) and (x,T) plane, respectively. The dashed line in the 
(A, T) plane represents the line of equilibrium second-order phase transitions on the critical curve and thus represents 
the A-line. This line terminates at the phase separation curve (solid line) at the critical end point E. CE is the line of 
the first-order phase transitions between the 3 He-rich and the 4 He-rich 'normal' fluids with the critical point C. At the 
point E the curve CE turns into the line of first-order phase transitions between the 3 He-rich 'normal' fluid and the 
4 Hc-rich superfluid. The phase boundary A(T) between the 3 He-rich and the 4 He-rich 'normal' fluid or the 4 He-rich 
superfluid (represented by a solid line) is expected to exhibit a singular curvature ~ \T — Tb|~ q as T approaches the 
end point temperature from above or below [24]. a is the critical exponent describing the specific heat singularity 
on the critical line below the point E. Since in mean-field theory a = 0, there is no nonanalyticity at the end point 
whithin the present approach. In the (x,T) plane (see Fig. 1(b)) the coexistence curve is smooth at T = Te on the 
3 Hc-rich side. 

For K/J < 2.01681 the instability point Pt lies outside the coexistence curve for the 3 He-rich 'normal' fluid and 
the 4 He-rich 'normal' fluid. Therefore, as T decreases below T t , the phase separation between the 3 He-rich 'normal' 
fluid and the 4 He-rich superfluid commences at the point P t on the critical curve; hence P t = A is a tricritical point. 
The phase diagram for K/J = 1.8 is shown in Fig. 2. In the (A, T) plane (Fig. 2(a)) the first-order transition line 
between the 3 He-rich 'normal' fluid and the 4 He-rich superfluid which starts at A terminates at a triple point D where 
it meets the first-order transition lines between the 3 He-rich 'normal' fluid and the 4 He-rich 'normal' fluid (curve CD) 
and between the 3 He-rich 'normal' fluid and the 4 He-rich superfluid. 

For even smaller value of the ratio K/J there are no longer two distinguishable disordered phases, i.e., the line 
DC in Fig. 2(a) has shrunk to zero. The critical point for coexistence between the 3 He-rich 'normal' fluid and the 
4 He-rich 'normal' fluid, which occurs at x c = 1/2 and T c = K /4, disappears. In Fig. 3 we present the phase diagram 
for K/J = 1. In the (A, T) plane it exhibits a very simple form (Fig. 3(a)). The A-line meets the first-order transition 
line between the 3 He-rich 'normal' fluid and the 4 He-rich superfluid at the tricritical point A. The lines meet with 
a common tangent, a feature characteristic of the mean-field approximation. In the (x,T) plane (Fig. 3(b)) at the 
tricritical point the critical line T s (x) has the same slope as the phase-separation curve on the 3 He-rich side. The 
emergence of the type of phase diagram shown in Fig. 3 from the one shown in Fig. 2 takes place at that value 
of K/J for which there is an equilibrium between the phase at the critical point C and an ordered phase yielding 
K/J w 1.4298. For K/J slightly less than this value K/J = 1.4298, the tricritical point is located at T t w 0.397 and 
Xt ~ 0.206. Whithin mean-field theory the 3 He-rich side of the coexistence curve has a plateau for 0.8 < x < 0.3, 
i.e., right below the tricritical point small changes in temperature lead to pronounced changes in the concentration 
of 3 Hc. As K/ J is reduced further the tricritical point shifts to larger values of x and smaller values of T. Also the 
shape of the coexistence curve changes; the plateau disappears and the concentration of 3 He inceases more uniformly 
with the temperature. 

A phase diagram like that of Fig. 3 resembles qualitatively the experimental one [14], for which one finds for the 
tricritical point Ta/T s (0) = 0.4 and xa = 0.669. In our model xa is, however, always smaller than 0.5. 




(3.22) 



whereas 




(3.23) 
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IV. MONTE CARLO SIMULATIONS 



For the Monte Carlo treatment of the model Hamiltonian given by Eq.(2.1) a 4 He atom is represented by the 
normalized spin vector 



for each lattice site i in the spirit of the standard XY model. A 3 He atom on lattice site i is represented as Si = (0, 0). 
Consequently the occupation number ti on lattice site i is given by ti = Sj| and the interaction cos(9i — 0j) between 
two 4 He atoms is given by the computationally more favorable scalar product cos(9i — 6j) = Si ■ Sj. The lattice is 
simple cubic with periodic boundary conditions and L lattice sites in each direction. The Monte Carlo algorithm is 
based on various standard procedures which we discuss briefly in the following. 

In order to explore the phase space of the model, two types of updates are needed: (i) spin flip updates and 

(ii) particle insertion and deletion updates. The spin flip updates are responsible for the creation of long - ranged 
magnetic order which in our model represents the normal - supcrfluid transition. This can be of first or second order 
depending on the concentration x = 1— < ti > of 3 Hc (ti = 0) in the system (see Fig. 3). The particle insertion and 
deletion updates are responsible for the demixing transition (phase separation) in our model. This transition can also 
be first or second order depending on the coupling constants in the model (see Figs. 1 and 2). In our simulation we 
are primarily interested in that regime of coupling constants, for which the phase diagram resembles that of actual 
3 He - 4 He mixtures (Fig. 3) and therefore the possibility of a second order (critical) demixing transition is not taken 
into account for the selection of the Monte Carlo moves. We therefore use the following methods in our Monte Carlo 
simulation: (i) single particle insertion and deletion and (ii) single spin flip according to the Metropolis algorithm [25] , 

(iii) single cluster spin flip according to the Wolff algorithm [26], and (iv) overrelaxation updates of the spin degrees 
of freedom at constant configurational energy [27]. For each particle insertion or single spin flip move the new spin 
state is randomly selected from the even distribution on the unit circle. The projection vector for the embedding part 
of the Wolff algorithm [26] is also chosen randomly from the even distribution on the unit circle. 

The above update methods are performed in sweeps over the whole lattice, where each spin flip sweep (ii), (iii), or 

(iv) is preceded by a Metropolis particle insertion and deletion sweep (i) . Cluster updates of the particle configuration 
according to the embedding algorithm of Ref. [26] are disregarded, because the critical demixing transition will not 
be explored here. 

The three basic Monte Carlo updates (ii) - (iv) outlined above are combined according to the hybrid Monte Carlo 
idea [28] to ensure efficient configuration space exploration also for second order (critical) transitions to long-ranged 
magnetic order. One hybrid Monte Carlo step consists of 10 individual steps each of which can be one of the updates 
listed above. The Metropolis and the Wolff algorithm work the standard way, in which the acceptance probability p 
of a proposed spin flip in the Metropolis algorithm is defined by the local heat bath rule 



where AE is the change in configurational energy of the proposed move. The overrelaxation part of the algorithm 
performs a microcanonical update of the configuration by sequentially reflecting each spin in the lattice at the direction 
of the local field, i.e., the sum of the nearest neighbor spins, such that its contribution to the energy of the whole 
configuration remains constant. The implementation of this update scheme is straightforward, because according to 
Eq.(2.1) the energy of a spin with respect to its neighborhood has the functional form of a scalar product. The form 
hybrid Monte Carlo step depends on the region of the phase diagram to be explored. In the vicinity of the first-order 
phase separation line typically six Metropolis (M), one single cluster Wolff (C), and three overrelaxation (O) updates 
are performed. The individual updates are mixed automatically in the program to generate the update sequence M 
M O M O M M O M C for the magnetic degrees of freedom. 

The random number generator is the shift register generator R1279 defined by the recursion relation X n = X n - p © 
X n - q for (p, q) = (1279, 1063). Generators like these are known to cause systematic errors in combination with the 
Wolff algorithm [29]. However, for lags (p,q) used here these errors are far smaller than typical statistical errors. 
They are further reduced by the hybrid nature of the algorithm [28] . 

The hybrid algorithm is well suited to explore second-order phase transitions. However, it is unable to overcome the 
exponential slowing down of the algorithms included in our hybrid scheme in the vicinity of a first-order transition, 
e.g., the first-order magnetic (i.e., superfluid) transition for higher concentrations of 3 He particles, i.e., for occupation 
numbers ti = in our model. In order to resolve this problem while keeping the benefits of the hybrid scheme we have 
embedded the hybrid Monte Carlo method in a simulated tempering environment [30]. According to the simulated 
tempering idea the temperature is treated as a random variable which performs a random walk inside a predefined 
temperature interval. In our simulation this temperature interval is represented by a discrete set of temperatures, 



Si = (cos 9i, sin 8i) 



(4.1) 



p(AE) = l/[cxp{AE/k B T) + 1], 



(4.2) 
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which are spaced closely enough to allow sufficient overlap of the corresponding energy distribution functions. The 
required rewcighting factors [30] are estimated from short runs, one for each pair of neighboring temperatures, and 
checked a posteriori by monitoring the probability distribution of the temperatures - which should be essentially flat 
- during the production run. Deviations of up to 20% from a flat temperature distribution are tolerated. 

The Monte Carlo scheme described above is employed for lattice sizes L between L — 12 and L = 60. For 
each choice of parameters we have performed at least 12 blocks of 10 3 hybrid steps for equilibration followed by 
another 10 3 hybrid steps to estimate the reweighting factors for each pair of neighboring temperatures and finally 
followed by 4 x 10 4 hybrid steps for measurements. The measurement block is controlled by an outer loop in which 
a new temperature is proposed according to the predetermined weight factors [30] after each hybrid Monte Carlo 
step. Apart from standard thermodynamic quantities the distribution functions of the total energy, the density, 
and the modulus of the magnetization arc monitored using histogram reweighting and extrapolation techniques [31] 
within the measurement block. Their statistical errors are estimated following standard procedures resulting from 
the statistical independence of different measurement blocks. Unless otherwise stated all error bars quoted in the 
following correspond to one standard deviation. They are displayed only when they exceed the symbol sizes. The 
simulations have been performed on DEC Alpha Workstations and Pentium III PCs. 

V. MONTE CARLO RESULTS 

Our primary interest in this study is to use Eq. (2.1) to model 3 He- 4 He mixtures in the tricritical region and we 
therefore restrict the numerical investigation of the statistical model described by Eq. (2.1) to the case J = K for 
which the phase diagram corresponding to this model Hamiltonian has the same topology as for the liquid phases of 
3 He- 4 Hc mixtures. The tricritical point marks the onset of demixing into a spin ( 4 He) rich fluid and a vacancy ( 3 Hc) 
rich fluid, where the spin rich fluid simultanously exhibits long-ranged magnetic order of the XY type (superfluidity). 

The phase diagram is most conveniently investigated by the inspection of distribution functions (histograms) for 
various thermodynamic quantities [32]. However, the computational expense of the method described in Ref. [32] 
in d = 2 is prohibitive in d = 3 for any appreciable system sizes. We therefore resort to a simpler though less 
accurate approach which allows us to treat larger systems and is accurate enough for our purposes. In the following 
all temperatures are measured in units of the critical temperature T s (0) of the XY model on a simple cubic lattice in 
d = 3, which is given by K c = J/(k B T s (0)) = 0.45415(5) [33]. The chemical potential A is measured in units of the 
magnetic coupling contsant J (see Eq. (2.1)). 

A. Order parameter distribution at tricriticality 

The key feature of the 3 Hc- 4 He phase diagram is the presence of a tricritical point. The task to locate the tricritical 
point for the model Hamiltonian given by Eq. (2.1) is aided by the observation that d = 3 is the upper critical 
dimension for tricriticality. It is therefore reasonable to assume that the distribution function of the magnetic order 
parameter essentially takes the mean-field (Landau) form. We will give some a posteriori evidence below that this 
assumption is indeed correct, but an accurate numerical proof of it is beyond the scope of this paper. 

The magnetic order parameter, i.e., the magnetization is defined by 

M = (M x ,M y ) = L~ 3 (5.1) 

i 

where U = 0, 1 characterizes the presence of 3 He or 4 He at lattice site i and Si = (cos0j,sin0j) is the standard spin 
variable of the XY model. In terms of the modulus m = |M| of the order parameter the distribution function P(m) 
is assumed to take the form 

P{m) = P m cxp(-^m 2 - Bm 4 - Cm 6 ) (5.2) 

according to Landau theory in the tricritical region, where the absence of symmetry breaking fields is assumed. The 
parameters A, B, and C essentially play the role of the Landau-Ginzburg model parameters (see Sec. VI below) and 
they depend on the temperature T and the chemical potential A (see Eq. (2.1)), where C is manifestly positive, but 
A and B may change sign. For system sizes L = 12, 18, 24, 36, 48 and 60 simulations have been performed along 
various paths in the (T, A) plane of the phase diagram and the data recorded for P(m) have been fitted according to 
Eq. (5.2) using Pq, A, B, and C as fit parameters. For each system size L a pseudo tricritical point (T t (L), A t (L)) has 
been identified by the requirement A = B = within the corresponding statistical error. It turns out, that Eq. (5.2) 
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indeed captures the shape of P(m) rather accurately over several orders of magnitude for P in the pseudo tricritical 
regime (see below). In particular, higher powers of m compatible with the symmetry such as m 8 can be safely ignored. 
Possible logarithmic corrections to P(m) [34] could not be identified from the numerical data unambiguously. We will 
comment on other logarithmic corrections later. 



B. The tricritical point 

From the procedure outlined above we obtain a sequence of pseudo tricritical points (T t (L), A t (L)) which can be 
extrapolated to the bulk limit L — > oo. In order to do this one has to identify the functional form of the L dependence 
of the pseudo tricritical point. Within our mean-field picture of the tricritical behavior of our model the coefficients 
A and B in Eq. (5.2) are given by the linear combination 

B \ ..( T-T t 
A )= M A-A t 



M a _ a. (5-3) 



in the vicinity of the tricritical point, where M. is the coefficient matrix. From Eq. (5.3) and finite-size scaling 
arguments one concludes that for sufficiently large L, T t (L) —T t and A t (L) — A t are governed by a linear combination 
of L~ dA and L~ ds , where d^ and d B are the scaling dimensions of the parameters A and B in Eq. (5.2) given by 
d,A = 2 and ds = 1 apart from logarithmic corrections [34]. We therefore arrive at the following functional form of 
T t (L) and A t {L): 

T t (L)=T t + | + ^ 

A t (L) = A t + 6 j- + ^, (5.4) 

where the coefficients t\, t 2 , 8\, and S 2 can be obtained from the inverse matrix and the finite-size relations 

A = aL- 2 and B = bL' 1 (5.5) 

for A and B evaluated at (T, A) = (T t (L) , A t (L)) for any system size L. The coefficients a and b are nonuniversal 
metric factors. Eq.(5.4) is used to fit the numerical data for T t (L) and A t (L) in order to obtain an estimate for the 
location of the tricritical point. Logarithmic corrections as given in Ref. [34] can be included in Eq. (5.4), but they 
are omitted here because the quality of fit does not change substantially when they are included. The results are 
shown in Figs. 4 and 5. The finite-size behavior of T t (L) and A t (L) is accurately captured by Eq. (5.4). Both the 
coefficients ti, t 2 and 61, 82 have the same sign and the second coefficient is substantially larger that the first one in 
both cases. Therefore both coefficients must be kept in order to obtain an acceptable fit. The quality of fit can be 
measured in terms of the reduced \ 2 which is 0.15 in Fig. 4 and 0.42 in Fig. 5. We thus obtain the extrapolated 
values 

T t /T s (0) = 0.7438(4), A t /J = 3.436(2) (5.6) 

as our estimate for the location of the tricritical point, where the statistical uncertainty affects the last digit by the 
amount given in parenthesis. In these units the coefficients in Eq. (5.4) are given by 

ti/T s (0) = 0.23 ± 0.02 , t 2 /T s (0) = 1.34 ± 0.24 , (5.7) 
81/ J = -1.05 ± 0.11 , 82/ J = -8.5 ± 1.4 . 

Another aspect of Eq. (5.4) is field mixing [32], because the finite-size corrections L _1 and L~ 2 are uniquely related to 
the coefficients (scaling fields) B and A in Eq. (5.2), respectively. In the vicinity of the tricritical point one therefore 
obtains from Eqs. (5.3), (5.4), and (5.5) by a matrix inversion 



( B/b\ = 1 / S 2 -t 2 \ f T - T t 
\ A/a ) tl 8 2 - t 2 <Si \Si h ) \ A - A 



(5.8) 



According to our mean field picture of tricriticality in d — 3 the coexistence line A = A(T) in the vicinity of T = T t 
should be associated with the line B = const. = in the vicinity of A = 0. If we linearize the coexistence line near 
the tricritical point according to 
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A(T) = A t + A; (T - T t ) 



(5.9) 



we obtain from Eqs. (5.7) and (5.8) for the slope A' t at the tricritical point 

A; = S 2 /t 2 = (6.4 ± 2.2) J/T a (0). (5.10) 

Despite its large statistical error this result serves as a valuable guideline for the further exploration of the phase 
diagram. 

The foundation of the above estimates is the quality of the fits of Eq. (5.2) to the measured order parameter 
distribution functions. We illustrate the quality of these fits in Fig. 6 for L = 36 at the corresponding pseudo 
tricritical point T = It (36) and A = A* (36). The shape of the distribution function P is essentially captured by Eq. 
(5.2) over more than three orders of magnitude. The parameters A and B vanish within their statistical errors. The 
reduced x 2 of the fit is 0.71. If A = B = is enforced, i.e, the fit is performed only with the parameters P and C, 
the reduced % 2 increases to 0.92. For all other system sizes investigated the situation is similar. We will return to the 
finite size behavior of P(m) after the discussion of finite-size scaling. 



C. Finite-size scaling 

A naive finite-size scaling ansatz for a thermodynamic quantity X(A, B, L) near a tricritical point in d = 3 is given 
by (compare Eq. (5.2)) 

X(A,B,L)^L- dx f x (AL 2 1 BL), (5.11) 

where fx(x,y) is the finite-size scaling function associated with the quantity X and dx is its scaling dimension. 
Logarithmic corrections have been disregarded for simplicity. For the sequence of the pseudo tricritical points (T, A) = 
(T t (L), A t (L)) one has A = aL~ 2 and B = bL^ 1 (see Eq. (5.5)). In this case one therefore expects X to display the 
scaling behavior 

X(aL- 2 ,bL-\L) = L- dx f x (a,b) = X L~ dx , (5.12) 

which can be conveniently checked numerically. However, near tricritical points in d = 3 one has to consider logarithmic 
corrections to the naive scaling and these have been examined in Ref. [34]. We therefore only quote the results 
corresponding to Eq. (5.12) for the average magnetization (to), the specific heat C and the magnetic susceptibility X. 
One obtains 

L\~ 1/2 / LX 1/4 



(m) = mo UJ [ ln r 



C = Co Tq [1. To ) , (5.13) 
L\ 2 / i^ 1/4 



x = x °[t ) IN 

where the nonuniversal amplitudes to , C , X and the length scale l are used as fit parameters. The corresponding 
results are summarized in Figs. 7-9. The data are compatible with the finite-size scaling behavior given by Eq. 
(5.13) (solid lines). The logarithmic corrections turn out to be essential. Disregarding these leads to pure mean field 
behavior which is not compatible with the data (dashed lines). For the specific heat displayed in Fig. 8 deviations 
from the expected behavior occur for larger system sizes L — 48 and L — 60. These may be due to the proximity 
to the first-order demixing transition, which is characterized by a finite latent heat. Including a finite background 
contribution to the specific heat as an additional fit parameter does not improve the fit. In particular, the attempt 
to fit pure mean field behavior to the data shown in Fig. 8 leads to a negative value for the background specific heat 
which is inconsistent with thermodynamics. The susceptibility shown in Fig. 9 appears to agree with the expectation 
for all system sizes whereas the average magnetization shown in Fig. 7 displays a deviation for L = 60. The values 
of l n obtained from the fits shown in Figs. 8 and 9 are consistent (l — 6.3 ± 0.5 and l — 6.2 ± 0.4, respectively) 
whereas lo = 1.3 ± 0.3 obtained from (to) according to Fig. 7 deviates strongly from the aforementioned estimates 
for Iq. One of the reasons may be that (to) depends rather weakly on L and l as compared to C and X . Therefore 
the estimation of lo from (to) is more susceptible to statistical or systematic errors in the magnetization data. The 
large relative error of the actual estimate l = 1.3 ± 0.3 seems to indicate this. Corrections to the leading asymptotic 
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behavior given by Eq. (5.13), which cannot be taken into account on our current data basis, may therefore also play 
a role. 

The scaling behavior of the order parameter distribution function P{m) within the scope of Eq. (5.2) is determined 
by the finite-size behavior of the parameter C. In order to compensate the finite-size effects induced by (to) wc define 
the effective coupling parameter 

C eff ^C(mf, (5.14) 

where C is taken from the fit of Eq. (5.2) to the distribution function data along the sequence of the pseudo tricritical 
points (A = B = 0) and (m) is taken from the fit of Eq. (5.13) to the magnetization data. The numerical result 
for C e ff according to Eq. (5.14) is displayed in Fig. 10 which shows a slow but systematic decrease of C e // with 
the system size. According to the renormalization group theory of tricritical behavior C e ff should play the role of 
the coupling constant at tricriticality, which is a dangerous irrelevant variable in d — 3 [34] . We therefore expect the 
finite-size behavior [34] 

C eff {L)=(c + c 1 ln^j (5.15) 

for the effective coupling parameter, where Iq = 6.2 is taken from Fig. 9 and Co and c\ serve as fit parameters to the 
data. The solid line in Fig. 10 displays this fit with Co = 9.9 ± 0.2 and c\ = 1.6 ± 0.1 and demonstrates, that the 
expected behavior according to Eq. (5.15) is consistent with the data. 

The degree of agreement between the finite-size scaling behavior observed and expected may also be considered 
as an a posteriori confirmation that the sequence of the pseudo tricritical points gives a reasonable estimate for the 
location of the tricritical point. However, some confirmation from a different source would still be desirable. 



D. Other distribution functions and the cumulant method 



In order to locate the first-order coexistence line one may inspect the distribution function P(n) of the particle 
density 

n = L- 3 ^2u. (5.16) 

i 

Near a first-order demixing transition P(n) displays two peaks, one at a higher density corresponding to the spin 
( 4 He) rich liquid and one at a lower density corresponding to the vacancy ( 3 He) rich liquid. At the tricritical point 
the two peaks coalesce and they separate increasingly as one follows the two-phase coexistence line towards lower 
temperatures. As a criterion to locate the coexistence line one may demand that the ratio of the statistical weights 
of the two liquids, i.e., the ratio of the areas under the respective peaks of P(n), should not depend on temperature. 
However, this criterion is only approximate, because a priori it is not clear what the value of the weight ratio should be. 
An accurate criterion can be obtained from the evaluation of the field mixing [32] (see Eq. (5.8)). From thermodynamic 
considerations one may determine a linear combination of particle and energy density such that the corresponding 
distribution function is symmetric on the coexistence line [32] . The generally unknown value of the weight ratio then 
has to be unity at coexistence. In principle, one may determine the correct mixing ratio of the densities from Eqs. 
(5.7) and (5.8). However, the statistical uncertainties of the coefficients given by Eq. (5.7) are too large for this 
purpose and their accurate evaluation is beyond the scope of this work. 

For the vector BEG model the task of locating the first-order coexistence line is aided by the observation that 
the spin rich fluid displays long-ranged XY type (superfluid) order when the demixing transition occurs. Apart from 
the particle density distribution (see Eq. (5.16)) we therefore also observe the distribution function of the magnetic 
energy density e m defined by 

s m = -JL~ 3 Mj cos(6 l - 9j). (5.17) 

<ij> 

For our choice J = K the demixing transition will also be indicated by a double peak structure of the distribution 
function P(e m ). Note that this will no longer be the case for sufficiently large K > J, for which the demixing 
transition precedes the onset of long-ranged magnetic order. By monitoring both distribution functions along various 
paths in the (T, A) plane of the phase diagram and by applying the constant weight ratio criterion to both we have 
redetermined the slope of the coexistence line in the vicinity of (T t , A t ) (see Eq. (5.6)) and found 
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a; = (5.o±o.i)j/r a (o). 



(5.18) 



Note that the new estimate given by Eq. (5.18) is consistent with the previous one given by Eq. (5.10). We 
furthermore observe, that the two peaks indeed merge into a single but broader one very close to the estimate of 
(T t ,A t ) given by Eq. (5.6). We illustrate this for P(n) in Fig. 11 for L = 36 along a straight path in the phase 
diagram according to Eq. (5.9) for the choice T t /Ts(0) = 0.7439, A t /J = 3.438 and A' t = 5.0J/T S (0) (see Eq. (5.18)) 
for three temperatures. In order to obtain a clear double peak structure in P(n) including the transition from and 
to a single peak along the chosen path a substantial amount of fine tuning for both T t and A t is required even for 
moderate system sizes. It is therefore very comforting that the values for T t and A t required to obtain Fig. 11 arc 
already within the error bars of the extrapolation estimate of the tricritical point given by Eq. (5.6). The structure 
of P(e m ) essentially looks the same so we do not reproduce it here. 

The location of the tricritical point along the coexistence line can be identified by a cumulant crossing of a suitably 
chosen density [32]. As we have not evaluated the field mixing here we use the cumulants of the magnetic energy 
density defined by Eq. (5.17) in order to investigate the cumulant crossing along the straight path used in Fig. 11. 
We define the Binder cumulant ratio for e m by 

j j _ -i (( g m ~ (gm)) 4 ) /r -i (\\ 

U <" = 1 - H{e m -(e m ) W (5 ' 19) 

The cumulant U £m as function of temperature for different system sizes is shown in Fig. 12. A unique crossing cannot 
be identified. However, the various crossings occur roughly where they are expected according to Eq. (5.6). If the 
smallest system L = 12 is excluded the crossings are located in the temperature interval 0.743 < T/T s (0) < 0.747 
which includes the estimate given by Eq. (5.6) near the lower bound. The crossings for larger systems tend to 
occur at lower temperatures. One of the reasons that a unique crossing does not occur is that both e m and n contain 
corrections to the order parameter of the demixing transition which can only be eliminated by a properly chosen linear 
combination of these quantities [32]. A second reason is insufficient fine tuning for larger systems, which becomes 
visible in the nonmonotonic behavior of U £m for L = 36 at lower temperatures which leads to a second intersection 
with U £m for L = 24. Despite the quantitative shortcomings of Fig. 12 the investigation of the cumulant crossing 
in combination with the other evidence presented above provides some independent confirmation that our initial 
assertion about the shape of the tricritical order parameter distribution function according to Eq. (5.2) is correct 
within the accuracy needed for the purpose of this work. 

Considerations in the spirit of Landau mean-field theory have played a major role for the analysis of our numerical 
data in the tricritical region. We therefore now turn to a detailed derivation of the Ginzburg - Landau model in the 
tricritical region of the vector BEG model. 



VI. LANDAU-GINZBURG MODEL FOR 3 HE- 4 HE MIXTURES. 

In this section we derive a two-parameter Landau-Ginzburg (LG) model describing bulk 3 He- 4 He mixtures near 
tricriticality. This derivation follows the construction of the <j> model for the standard critical phenomena from the 
Ising model. 



A. Derivation of the model 



Our starting point is the modified VBEG model for which instead of continuum orientations of the spin vector 
Si = (cos#j, sin#i) we consider L discrete orientations 6*p — 2irl/L, I — 1, uniformly distributed over the 

unit circle with L — > oo. With each orientation sf^ = (cos9f\ sin 9^) at the site i we associate the density tij with 
discrete values or 1. The total density of 4 He at the site i is 

L 

U = J2ki- (6-1) 

i=i 

As in the VBEG model we consider the close-packing case in which a 3 He atom at site i corresponds to tj = and 
a 4 He atom to t$ = 1. Thus a lattice site i is either occupied by a 3 He atom (tj = 0) or a 4 He atom associated with 
one of the L orientatons (U t i = 1 for I = lo, = otherwise, so that U = k,i — !)■ The Hamiltonian of this 

effectively (L + 1)- component mixture has the form 



12 



L L 

n=-J2i J J2J2 hitjM l) ■ Sp + Ktitj} + (6.2) 

<ij> 1=1 l' = l i 

where sj° • sf ) = cos^f - dp). 

In the Landau- Ginzburg model the effective Hamiltonian, depending on the local order-parameter fields, is obtained 
as a result of coarse-graining procedures. The procedure which gives an exact functional representation for the 
partition function for the corresponding microscopic Hamiltonian is the Hubbard-Stratonovitch transformation. The 
application of this method, however, is limited to microscopic Hamiltonians that can be expressed as a quadratic 
form. Here we use another approach [35]. 

Within the standard mean-field treatment of the lattice gas mixture defined by the Hamiltonian in Eq.(6.2) the 
ensemble-avaraged occupancy pij —< t it i > of the site i is obtained by minimizing the grand canonical function 

i 

at fixed T and A. The resulting minimum of Vt MF equals the equilibrium grand potential Vt$. The ideal or noninter- 
acting free-energy density for a (L + 1) component mixture on the close-packed lattice is [18] 

L 

fid(pi,l) = k B T{(l - p.) ln(l - p t ) + ^2pu In /),,;}, 

i=i 

with p l = Ya=i Pi,l- 

In the following we shall treat pi and p,,; as coarse-grained order parameter fields and adopt the mean field grand 
canonical function Vi MF form for the free energy of a particular local configuration of the order parameter. In the 
spatially uniform and orientationally disordered phase the equilibrium values of pi y i are constant and denoted by 
< Pi,i >— Q/L so that < pi >— Q. The actual values fluctuate around these mean values: 

Pi = Q + <t>t (6.5) 



(6.3) 



(6.4) 



Pi,i = j + f + Api,i = ^+A Piil , (6.6) 

which implies < <f>i >— and Apij — even without taking the thermal average. The fluctuation of the density 
Pi t i at the site i consists of an orientationally uniform part cf>i/L related to the fluctuation <\>i of the total He density, 

which is the same for all orientations, and a contribution Ap^; as an excess density of 4 He in the orientation S^'. 

Assuming Ap iy i and 4>i to be small we expand Q MF (pi, p iy i) in power series of the fluctuation fields Ap iy i and 4>i 
about the equilibrium value £l^f F (pi = Q, P ij = Q/L). Since we aim for deriving an effective Hamiltonian describing 
bulk 3 He- 4 He mixtures near tricriticality, in the expansion we keep terms to the sixth order in Apij and to quadratic 
order in A standard Taylor expansion gives 

6 

n MF ( PuP u) - n* iF (Q, q/l) = ^f F {4>^Ap^) 

The contribution linear in the fluctuation fields Q.^ IF vanishes since we expand f2 Mi * around its minimum. The other 
terms are: 

L L 

n^ F ( ( f >u Ap hl ) = -^2 {KMj + jY,Y, &Pi,i&PoA l) ■ spy 

<i,j> 1=1 l' = l 



(6.7) 



(6.8) 



^(^,A a 0--^EE{^( A ^) 2 + ^(Apu) 3 }, (6.9) 
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nffa, A Piil ) = *^2g{|^(A ftil ) 8 + ^^(Apu) 2 + ^(W>, (6.10) 
nr Aft,,) = ~Y. £{^<M W + ^ W + ^ (^,i) 6 }, (6-11) 

i (=1 

«r Aft,,) = ^££{^ ( Ap M ) 5 + i£#(W + ^7s(W}- (6.12) 

i 1=1 

The excess density in the orientation is a periodic function of Z with period L. Therefore, it can be expanded 
into a discrete Fourier series 

^4&' ieWil " (6 - i3) 
fc=i 

where we have chosen (Q/2L) as a normalization constant. The term corresponding to k = is excluded from the 
expansion due to the constrain Y^i=i A/Oj,; = 0. Since Apij is a real function, the Fourier components u it k and u^l-h 
are related by u* fc = u^L-k- Now, we neglect higher modes in the Fourier expansion (6.13) and approximate the 

excess density in the orientation sf^ by 

Api,i « ^ (u i;ie ^/^ + Ui^-ie-*^)') . (6.14) 
Using u* k ~ Ui^L-k and expressing tt^i in terms of its amplitude \ui\ and a phase Q we write Eq. (6.14) as 

Apm - (e'W^+He-'tW 1 "^') = f H cos((2,r/L)Z + 6) - fsf • u,, (6.15) 

where = \m\ exp(i£j). 

The approximation (6.15) is a 'coarse-graining' procedure which reduces degrees of freedom. The L — l independent 
quantities describing the oricntational degrees of freedom are replaced by a 2-component vector field u^. 
We define the effective Hamiltonian for the order parameter fields <j> and u as: 

n e " ui) = n MF ( Pi , pij) - n™ F {Q, q/l). (e.ie) 

Using Eq. (6.15) in Eqs. (6.8-6.12) we obtain 



n*ff = agt + n<l% (6.17) 



with 



<j,i> (=i i'=i 

L 



2 



L 



n 



til = ^ £{-4 E( g f 5 • u *) 2 + ^ E( s < (0 • + £( s f } • u *) 4 

i 1=1 ^ 1=1 1=1 

- • u *) 4 + # D s i° • u *) 4 + if B s " } • u *) 6 }- ( 6 - 19 ) 
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We note that in the above expression terms containing sums over all orientations of odd powers of (s|^ • Uj) vanish. 

As the next step we take the limit L — > oo. This ammounts to replacing j- ^2f =1 by J'q 2 ^ d6 and leads to the 
following relations (sf^ = (cos 6f\ sin Of^)): 



§£(S« -u,) 2 - £ jfVW -u,) 2 = §|u,| 2 , (6.20) 



|p S f) -a,) 4 - I^VW -u,) 4 = f |u,| 4 , (6.21) 
fX> ( ° - |^^f(sf) . Ui )« = §|u,| 6 , (6.22) 



and 



? 2 ^ 

i"2 



^EE^^sf-u^'-sf 



;=i z'=i 

Q2 ,2* ... ,2. 

4 



r^f) / W ^f) ( sf). Ui )(sf .u,)sf -Sf 
o Jo 

Ui ■ uj. (6.23) 



Finally, we asume that the fluctuating fields (f>i and vary slowly on the length scale of the lattice constant a. The 
continuum limit is obtained by considering a — > 0, considering z as a continuous variable r, and fa and turning 
into far) and u(r), respectively, while keeping the total volume V — a 3 N fixed. In this limit, one has 

-> / (6-24) 

For / being the smooth continuation to continuous arguments of a function defined on a lattice we use the following 
approximations (a — ► 0): 

^ ^ c\ £ 2 d o2 



and 



A;=l fe=l fe=l K 



1 d 2 f 



£{/(r + ae k ) + /(r - ae fe )} - 2rf/(r) + a 2 ]T + • • • , (6.26) 

k=i k=i 



where efe , k = l,...,d are the unit lattice vectors. Thus 

d 



(6.27) 



E fifi - I dh E /( r M/( r + ae fe) + /( r - ae *)} 

<i,j> " ^ fe=l 

_^I a -3 y d 3 r{2d/ 2 (r) - a 2 (V/) 2 }. 
As a result Eq. (6.16) is replaced by 

= tf[n G + n int ] (6.28) 
with the Gaussian contribution flc, in which the fields <f> and u are uncoupled, 
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G = J dv{ l -a 1( f + i(V^) 2 + ia 2 |u| 2 + ic(Vu) 2 }, (6.29) 

and the interaction contribution, which couples <f> and u, 

flint = J dr{r^\u\ 2 + a 12( p 2 \u\ 2 + b\u\ 4 + r 2 0|u| 4 + 6 12 2 |u| 4 + e|u| 6 }, (6.30) 

where we have chosen the length unit such that a = 1. 

The coupling constants in the effective functional are given explicitly by: 

k B T 1 

fll = -?TQ(l— Q)- Z > (6 ' 31) 
fc B TQ zQ 2 J 

tt2 = T2-T K> (6 ' 32) 
O 2 J 

1 k B T 3 fc B T 

ri = -4^T< r2 = -32^T' (6 ' 34) 

6=«*^, e JM, ( 6 .35) 

32 K ' 96 X ' V ; 

fll2= 4Q-?T' 6l2 = T6Q-?T- (6 ' 36) 

Equations (6.29), (6.30) and (6.31-6.36) define the Landau- Ginzburg model for 3 He- 4 He mixtures in terms of ther- 
modynamical quantities and two parameters, J and K, characterizing the system. 

B. A-line and tricritical point 

In this subsection we determine the A-line and a tricritical point within mean-field theory for the LG model derived 
in the preceding subsection. To this end we consider spatially uniform order parameter fields. 

Mean-field theory amounts to approximating the thermodynamic free energy by the minimum of the effective 
Hamiltonian which corresponds to the saddle point path contribution to the partition function: 

I3F MF = min 0iU /3O e ^[0, u]. (6.37) 

The mean-field solution for 0(r) is determined by 

fk= a ( " 8) 

For spatially uniform fields the above minimum condition yields the following relation between <j> and |u|: 

0(oi + 2a 12 |u| 2 + 26 12 |u| 4 ) + n|u| 2 + r 2 |u| 4 = 0. (6.39) 

Near a tricritical point both the fields <j> and |u| are small. Therefore it is sufficient to consider only a linear coupling 
between |u| 2 and <j>, neglecting the higher order terms: 

^-^|u| 2 + 2^|u| 4 + 0(|u| 6 ). (6.40) 
Inserting Eq. (6.40) into Eqs. (6.29) and (6.30) we obtain 
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^-^ 2 |u| 2 + (6-^)|u| 4 + e '|u| 6 , (6.41) 

with 

2 

e' = a l2 r -\- — +e. (6.42) 

The condition a 2 — yields the equation for the critical line which is in agreement with Eq. (3.17). If a\ is negative, 
e' is positive and there is a tricritical point determined by 

1 r 2 

a 2 =0 b--— = 0. (6.43) 

2 CL\ 

The solution of these two equations coincides with the expressions given by Eqs. (3.19) and (3.20), i.e., the tricritical 
point of this LG model is located at the same temperature and concentration of 3 He atoms as the tricritical point in 
the VBEG model studied in Sec. Ill within mean field approximation. 



VII. SUMMARY 



By using molecular-field approximations and Monte Carlo simulations we have investigated a three-dimensional 
version of the generalized spin-1 Mume-i?mery- Griffith model (Eq. (2.1)) of 3 Hc- 4 Hc mixtures with a two-component 
vector order parameter, mimicing the phase of the wavefunction of 4 He atoms. This work is a first step to study the 
Casimir force and other surface and finite-size effects in 3 He- 4 He mixtures films near their tricritical point. We have 
obtained the following main results: 

1. The topology of the phase diagram depends on the ratio of the interaction parameters K/J, where K is related 
to the "He-^He interactions (Eq. (2.2)) and J to the superfluid density (Eq. (2.4)). There are three different 
types of the phase diagram, which are similar to those found in the BEG model within the molecular-field 
approximation. For large values of K/J, i.e., for K/J > 2.01681, there exist three different phases: a 3 He-rich 
'normal' fluid, a 4 He-rich 'normal' fluid, and a 4 He-rich superfluid (see Fig. 1). As the temperature is lowered, 
the mixed normal fluid phase separates into two 'normal' fluids differing by the concentration x of 3 He. This 
phase separation ends at a critical point. At lower temperature, the phase separation is into a superfluid and 
a 'normal' fluid. The A- line T s (x) of second-order phase transitions between a 4 He-rich 'normal' fluid and a 
4 Hc-rich superfluid terminates at the 4 He-rich branch of the phase-separation curve at the critical end point. 
This 'critical end-point' type of the phase diagram was the only one found in previous studies [15,16] of the two- 
dimensional version of the model. In three dimensions we find two additional topologies of the phase diagram 
as the ratio K/J is decreased. For 1.4298 < K/J < 2.01681 the phase diagram is the richest (see Fig. 2). As in 
the previous case, there are three different phases and a critical point of the phase-separation curve between two 
'normal' fluids differing by the concentration of 3 He. In addition, there is a tricritical point at the end of the 
A-line beyond which a first-order phase transition between a 4 He-rich superfluid and a 4 He-rich 'normal' fluid 
takes place. There is also a triple point at which three different phases coexist at different concentrations. For 
K/J < 1.4298 the phase diagram simplifies (see Fig. 3). There is no longer a 4 He-rich 'normal' fluid phase and a 
critical point. The A-line meets the first-order phase separation line between 4 He-rich superfluid and a 3 He-rich 
'normal' fluid at the tricritical point. The A-line is given by T s (x) = zJ(l — x)/2, where z is the coordination 
number of the lattice. The temperature of the tricritical point is T A /T s (0) = (1 + 2K/ J)/ (2 + 2K/ J) (Eq. (3.19)) 
and the concentration xa of 3 He at this point is given by Ta/T s (0) = 1 — xa (Eq. (3.20)). This type of the 
phase diagram is similar to that observed experimentally, although in our model xa is always smaller than 0.5 
whereas x e A xp = 0.669. 

2. The existence of the tricritical point is confirmed by Monte Carlo simulations and in the units (T t /T s (0), A t / J) it 
coincides with the mean-field prediction remarkably well (see Fig. 3). At the tricritical point the order parameter 
distribution function takes its mean-field form, where the presence of logarithmic corrections could not be 
excluded within the accuracy of the existing data. On the other hand finite-size scaling of several thermodynamic 
quantities reveals the presence of logarithmic corrections in accordance with theoretical expectations. The two- 
phase coexistence line in the (T, A) plane of the phase diagram has been determined form a constant weight 
ratio criterion for energy and density histograms. The location of the tricritical point is also indicated by a 
crossing of the cumulant ratio of the magnetic portion of the energy density measured along the coexistence line 
(see Fig. 12). We conclude that mean-field theory provides a reliable approach for studying the VBEG model 
in d = 3. 
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3. Starting from the VBEG model and discretizing the orientations of the spin vector we have derived the continuum 
Landau-Ginzburg model for 3 He- 4 He mixtures near the tricritical point encompassing the concentration field <f> 
and a two-component vector field u corresponding to the orientational order. In the effective Hamiltonian we 
consider the modulus of u up to its sixth power and the field 4> U P to quadratic terms, which is sufficient to 
study a tricritical point. The coupling constants appearing in this Landau-Ginzburg theory are given explicitly in 
terms of thcrmodynamical quantities, the temperature, the mean concentration of 4 He, and the two interaction 
parameters J and K characterizing the VBEG model. Mean-field theory for this LG model yields the same 
equation for the critical A-line as the molecular-field approximation for the lattice VBEG model. The LG model 
provides a linear coupling between |u| 2 and <p which yields the same coordinates of the tricritical point as the 
lattice VBEG model. 
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FIG. 1. Phase diagram in the (A,T), (a), and (x,T), (b), plane for the model given by Eq. (2.1) obtained whithin mean-field 
theory for K/J — 2.4. x is the 3 He concentration. There are three phases which can be identified as a 3 He-rich normal fluid 
(M — 0,x = 1 — Q large), a 4 He-rich normal fluid (M = 0, x small), and a 4 He-rich superfluid S (M 0,x small). In (a) 
the dashed line represents second-order phase transitions and corresponds to the A-line; full lines are the loci of first-order 
phase transitions. For this value of K/J there is no tricritical point. The A-line of second-order phase transitions termi- 
nates at the phase-separation curve at the critical end point E. The two-phase region in (b) and the line of first-order phase 
transitions in (a) end at a critical point C. The coordinates of the critical points are: C = (A/T s (0) = 2.4, T/T s (0) = 1.2), 
E = (A/T s (0) = 2.4, T/T s (0) ~ 0.89) and C = (T/T s (0) = 1.2, x = 1/2) and E = (T/T s (0) w 0.89, x « 0.107). In (b) the A-linc 
is given by T a (x)/T a (0) = 1 — x (see Eq.(3.17)). The two-phase region in (b) between C and E is symmetric about x = 1/2. 
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FIG. 2. Same as in Fig.l for K/J — 1.8. For this value of K/J the A-line ends at a tricritical point A beyond which there 
is a first-order phase transition between the 4 He-rich superfluid S and the 4 He-rich normal fluid. At even lower temperatures 
there is a triple point D. The coordinates of these points are: A = (A/T s (0) ~ 1.704, T/T s (0) « 0.821), C = (A/T s (0) = 1.8, 
T/T s (0) = 0.9), D = (A/T s (0) = 1.8, T/T 3 (Q) « 0.793) and A = (T/T s (0) ~ 0.821, a; « 0.179), C = (T/T s (0) = 0.9, a; = 1/2) 
and D = (T/T s (0) = 0.793, xi w 0.154, £2 w 0.216, X3 w 0.784). In (b) there are two-phase coexistence regions below A and 
below C which join for three-phase coexistence at D (dotted curve). 
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FIG. 3. Same as in Figs. 1 and 2 for K/J — 1. The A-line T s (x) and the first-order phase separation line meet at the 
tricritical point A. The 3 He-rich 'normal' fluid phase is denoted by N. In (b) Monte-Carlo data for the phase boundaries are 
indicated by pluses which are connected by thin lines representing the Monte Carlo phase boundaries. The inset shows the 
results on an expanded scale near the tricritical point A. The coordinates of the tricritical point within mean-field theory are: 
A = (A/T s (0) « 0.776, T/T s (0) = 0.75) and A = (T/T s (0) = 0.75, £ = 0.25). In (b) the tricritical point A as obtained from 
Monte Carlo data is also denoted by a dot and has the coordinates (T/T s (0) = 0.744, x = 0.26). 
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FIG. 4. Pseudo tricritical temperature T t (L) (x) vs. 1/L measured in units of the critical temperature T s (0) of the XY 
model on a s.c. lattice in d — 3 for J — K (see Eq. (2.1)). Error bars correspond to one standard deviation. The solid line 
shows the fit of Eq. (5.4) to the numerical data. The arrow indicates the extrapolated value T t (see main text). The reduced 
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FIG. 5. Pseudo tricritical chemical potential A t (L) (x) vs. 1/L measured in units of the coupling constant J (see Eq. (2.1)) 
for J = K. Error bars correspond to one standard deviation. The solid line shows the fit of Eq. (5.4) to the numerical data. 
The arrow indicates the extrapolated value At (see main text). The reduced \ 2 of the fit is 0.42. 



27 




28 



0.50 




FIG. 7. Least square fit of Eq. (5.13) (solid line) to the simulation data for (m) (x). A fit to pure mean field behavior is 
shown for comparison (dashed line). Data and fit are normalized to the amplitude mo and lo = 1.3 ± 0.3. For L — 18 and 
L = 60 the data points deviate from the fit curve (solid line) by an amount larger than the statistical error. 



29 




FIG. 8. Least square fit of Eq. (5.13) (solid line) to the simulation data for the specific heat C (x). A fit to pure mean field 
behavior is shown for comparison (dashed line). Inclusion of a background contribution to C as an additional fit parameter 
does not improve the fit. Data and fit are normalized to the amplitude Co and l — 6.3 ±0.5. The deviations from the expected 
behavior (solid line) for larger systems may be due to the vicinity of the first-order demixing transition. 
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FIG. 9. Least square fit of Eq. (5.13) (solid line) to the simulation data for magnetic susceptibility X (x). A fit to pure 
mean field behavior is shown for comparison (dashed line). Data and fit are normalized to the amplitude Xq and l = 6.2 ±0.4. 




31 



FIG. 10. Effective coupling parameter C e // according to Eq. (5.14) (x) and a least square fit of Eq. (5.15) to the data 
(solid line). The observed decrease of C e ff is compatible with the logarithmic behavior of a dangerous irrelevant variable at 
the upper critical dimension. The reduced \ 2 of the fit is 0.16. 




FIG. 11. Particle density distribution P(n) for three temperatures along a straight path given by Eq. (5.9) in the tri- 
critical region as proposed by Eq. (5.6). The temperatures chosen are T/T s (0) = 0.7479 (x), T/T s (0) = 0.7439 (+), and 
T/T s (0) = 0.7399 (*). The parameters of Eq. (5.9) are T t /T s (0) = 0.7439, A t /J = 3.438, and A' t = 5.0J/T S (0). 
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FIG. 12. Cumulant ratio U Em according to Eq. (5.19) as function of temperature along the straight path used in Fig. 11 
for L — 12 (x), L — 18 (+), L = 24 (*), and L = 36 (□). Pairs of symbols are connected linearly to guide the eye. A unique 
crossing cannot be identified (see main text). 
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